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Abstract 

We discuss entanglement generation in a closed system of one or two atomic quantum dots 
(qubits) coupled via Raman transitions to a pool of cold interacting bosons. The system exhibits 
rich entanglement dynamics, which we analyze in detail in an exact quantum mechanical treatment 
of the problem. The bipartite setup of only one atomic quantum dot coupled to a pool of bosons 
turns out to be equivalent to two qubits which easily get entangled being initially in a product state. 
We show that both the number of bosons in the pool and the boson-boson interaction crucially 
affect the entanglement characteristics of the system. The tripartite system of two atomic quantum 
dots and a pool of bosons reduces to a qubit-qutrit-qubit realization. We consider entanglement 
possibilities of the pure system as well as of reduced ones by tracing out one of the constituents, 
and show how the entanglement can be controlled by varying system parameters. We demonstrate 
that the qutrit, as expected, plays a leading role in entangling of the two qubits and the maximum 
entanglement depends in a nontrivial way on the pool characteristics. 
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I. INTRODUCTION 



Entanglement is a hallmark of quantum mechanics and is fundamentally different from 
any correlation known in classical physics. It has become indispensable in quantum computa- 
tion because of its enormous capabilities to process information in novel ways l| . Entangled 
states form a base for quantum communication protocols such as superdense coding [2| or 
quantum teleportation 3|. 

Entanglement also plays a fundamental role in the physics of condensed matter and 
strongly correlated materials j4|. BCS wave function jsl, Laughlin ansatz j^, Kondo singlet 
01 are examples of highly entangled ground states. Entanglement entropy provides signifi- 
cant iiisight into quantum critical phenomena sj-l^ and is used to characterize topological 



order flli 



12], which is not described by the standard Ginzburg-Landau theory. 
It would be interesting for both quantum information and condensed matter if one could 
generate particular entangled states in a controlled manner. Cold atoms in optical lattices 
have become one of the favorite systems to tailor various many body phenomena, including 
multi-particle entanglement 13 1, which may be used for quantum computing [l^. Schemes 
based on arrays of cold atoms have also been proposed for generating highly entangled cluster 



states and entangling gates between distant qubits [ISHlTl]. The main advantage of cold 
atomic systems is that they are clean, well controlled and almost dissipationless. Moreover, 
impressive recent advances in fluorescence imaging have made it possible to probe with high 
idelity the on-site number statistics of cold atoms in optical lattices in the Mott regime 
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ll9|. 



In this paper we consider the emergence of entanglement in a system of two atomic 
quantum dots (AQDs), which constitute two qubits, coupled by optical transitions to a pool 
of bosons containing a finite number of particles. The bosons in the pool are interacting, and 
we consider the ratio of their interaction energy U with respect to the coupling T between 
the AQDs and the bosons as a main parameter in our problem. We first consider an even 
simpler system of just one AQD coupled to the pool and show that such a set-up is in fact 
equivalent to a two-qubit system whose ground state and excited state are a singlet and 
triplet state, respectively. We observe dynamic formation of entanglement in the system, 
starting initially from a product state. We show that the state can become maximally 
entangled in the regime U/T <^ 1. 
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We then proceed to the case of two AQDs coupled to the pool. Since the system is closed 
and the number of particles is conserved, this is equivalent to a qubit-qutrit-qubit realization, 
which is interesting because qutrits have larger entanglement capacitance than qubits 20 |. 
Although there is no well-defined entanglement measure for a tripartite system, we calculate 
a concurrence which tells us whether or not the state is entangled. We then trace out the 
qutrit (pool of bosons) and calculate entanglement of formation 2l| of the bipartite system. 
We demonstrate that dynamics of the entanglement of formation is consistent with that of 
the concurrence and discuss best entangling possibilities depending on U/T. Finally, we 



trace out one of the qubits and calculate both entanglement of formation and negativity 
for a remaining qubit-qutrit state. 



22| 



Our work is further motivated by a previous result by some of us [23| , in which we found 
that strong correlations can emerge between the two qubits coupled to a Bose-Einstein 
condensate. However, we have not found any entanglement in that case, most likely since 
the Bose-Einstein condensate was described by the Gross-Pitaevskii wave-function, which 
is effectively a classical field. The question we want to address here is therefore, whether 
quantum entanglement emerges if the BEG is not in the Gross-Pitaevski regime. 



II. MODEL: TWO AQDS COUPLED TO A POOL OF INTERACTING BOSONS 

Our set-up is displayed in Fig. 1. We consider a system of two atomic quantum dots 
coupled to a pool of bosons trapped in an external potential V^os- The system Hamiltonian 
reads 

H = Hhos + Haqd + H^ix, (1) 

The energy levels of Vbos are widely separated, so that particles can only occupy the lowest- 
lying level. This is effectively a one-site optical lattice in the Mott regime, so that iJ^os 
reads 

Hbos = Ei,n + ^n{n-l), (2) 

where Ei, is the energy of the bosons, n is the number operator and U is the interaction 
energy between the atoms. 

Each atomic quantum dot constitutes a qubit and can be created in the following way: An 
atom of a different hyperfine species from the bosons in Vfeos is trapped in a tight potential 
Vi^2- The interaction between these atoms is assumed to be so large that double occupancy 
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FIG. 1: Two atomic qubits coupled to a finite number of cold interacting bosons. The qubits are 
trapped by the tight potentials Vi and V2, the bosons are confined in an external potential Vhos- T 
is the optical transition assisted coupling between the qubits and the bosons. 



of the lowest lying level of Vi^2 is forbidden. The coup 



the pool is produced in an optical way suggested in |25|: An external blue-detuned laser 



ing between qubit and bosons in 



induces Raman transitions between the dot and the atoms in the pool. The possibility of 



such a scheme was also confirmed in a numeric simulation 26[. Importantly, the magnitude 
of the coupling between different atomic states can be rather large, since it is controlled by 
adjustable Rabi frequency of the external laser. However, the overlap between the wave- 
functions of the bosons in the trap Vbos and bosons of atomic quantum dots is considered to 
be small, so that density-density interaction between particles can be neglected. 

Under above-mentioned conditions we can map the bosonic creation (annihilation) op- 
erators (b) of the dot onto their Pauli matrix equivalents: 6^ — a+ and 6 — )■ a"_. The 
corresponding Hamiltonian is then 

2 
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^AQD = -|^$^(l + a,(z)), (3) 



i=l 

The argument of the Pauli matrix i is referring to the left i = 1 or the right i = 2 qubit. 
Here, A is the detuning energy necessary to avoid spontaneous emissions during the Raman 
process. 

Finally, the interaction between the qubits and interacting bosons in the pool reads 

H^i, = T [aa+{l) + a^a^l) + aa+{2) + ata_(2)] , (4) 
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where (a) creates (annihilates) a boson in the pool and the coupling T is proportional 
to the Rabi frequency VLr of the above-mentioned Raman transition, i.e., T ~ hVlR. We 
consider a closed system and the total particle number N = n + \{1 + 0-2(1)) + |(1 + o-z{2)) 
is conserved, i.e., [H,N] = 0. In the following we tackle the problem of the dynamics of 
entanglement in this tripartite system. 



III. A SINGLE AQD COUPLED TO THE POOL: BIPARTITE ENTANGLEMENT 

It is instructive to start from an even simpler bipartite case of only one AQD coupled 
to the pool of bosons. Due to the particle number conservation the corresponding Hilbert 
space "H = "Hbos ® 'Haqd is spanned just by two orthonormal states 

|n,0) = \n) ® |0), 
|n-l,l) = |n-l)®|l). (5) 

The state |0) corresponds to "spin-down" state of a dot, while the state |1) corresponds to 
a "spin-up" . The wave function in this case is 

\^{t)) = c,{t)\n,0)+C2{t)\n-l,l) (6) 

with the normalized coefficients |ci(t)p -|- |c2(t)p = 1. The wave function ([6]) is in fact 
equivalent to a standard wave function of two spins- 1/2 or two qubits, with two of the 
four coefficients equal to zero from the onset. We can immediately write down the time- 



dependent concurrence C(t) of 
as an entanglement measure 



_;he state as C{t) = 2|ci(t)c2(t)| and use it in the following 
21| . The state is entangled whenever both coefficients ci and 
C2 are nonzero. Note that the initial state |^I^(0)) = |n,0) is always unentangled. 
The Hamiltonian in the basis ([5]) reads 

H, = H[h+\ ^ , (7) 

-E-U{n-l) j 

where H[ = Ef,n + ^Un{n — 1) is the diagonal part of the Hamiltonian (which we neglect in 
the further consideration) and E = + A is the energy difference between the traps. The 
Hamiltonian (JTj) describes a two-level system with Rabi frequency ur = T^fn and detuning 



b = ^[—E — U{n — 1)], so that the eigenenergies are E± = 6 ± + ujji- We note that for 
6 = the ground state is a singlet and the excited state is a triplet. 



We now consider the time evolution of the concurrence C{t), or equivalently, the entan- 
glement of the state ([6]), initially unentangled |^(0)) = \n,0). First, we discuss the case of 
small detuning 6/T ^ 1. The maximum concurrence is then always close to unity and the 
concurrence is approximately described by C ~ \exp{—2i6t) sm{2^/nt)\, shown in Fig. 2(a) 
for various U/T and fixed E and n. In order to analyze the results, we introduce a con- 
venient quantity tent referred to as "entanglement period", which is the time span between 
two subsequent instances of zero entanglement. In Fig. 2(a), for example, tent = 0.5 for 
U/T = and t^nt = 1 for any finite U /T (time is in the units of 1/T). 

In Fig. 2 {b) we plot the entanglement period and the maximum concurrence Cmax as 
a function of the interaction U/T for two values of n. We see from the results in Fig. 2 
(b) in order to increase the entanglement period and maximally possible entanglement one 
should either decrease U /T or the number of particles in the bosonic pool n. Large values 
of n as well as large values of U/T do not generally advocate entanglement in the system. 
This happens because the amplitude of the Rabi oscillations scales as ur/ {uj^ + 6"^), so that 
large values of 6 contribute to the suppression of entanglement. 



IV. TWO AQDS COUPLED TO THE POOL: BIPARTITE AND TRIPARTITE 
ENTANGLEMENT 



We consider now our initial tripartite system: two AQDs coupled to the pool of interacting 
bosons. The corresponding Hilbert space is a tensor product of three subspaces "H = Hbos ® 
'HaqDi ®'Haqd2 a-iid is spanned only by four basis vectors (because of particle conservation) 
which we choose in the following way 



{|n + 1,00), |n,01), |n, 10), |n- 1,11)}. 



(8) 



Note that the pool is effectively a three state system, i.e., a qutrit, described by the states 
"spin-up" |?7, + 1), "spin-fiat" \n) and "spin-down" \n — l). 

We can now write the Hamiltonian of the tripartite system in the basis ([8]) as 

I E + Un T^/^T+l T^^^^Tl \ 

TV^n^ 

TV^J^TT Ty/n 

\ T^ -E + U{l-n) J 



H = Hnl 



0-'4 



(9) 
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FIG. 2: (a) Time evolution of the concurrence of the bipartite setup for n = 10, E/T = 0.01 and 
the interaction U /T = (black), U /T = 0.06 (red) and U /T = 0.1 (blue), {h) Dependence of 
the entanglement period ignt (red curves) and maximum concurrence Cmax (blue curves) on the 
interaction U /T for n=10 (circles) and n=4 (diamonds) bosons in the pool. 

where Hq = Etn + \Un{n — 1) — A is the diagonal part. One of the eigenstates of this 
Hamiltonian is a singlet in the AQD sector |^s) = ® (|01) — |10)) corresponding to 

an eigenvalue equal to zero. All other eigenvalues have to be determined numerically from 
the matrix ([9]). We just mention that the singlet is the second excited state, which is very 
close to the third excited state (the closer the smaller the interaction energy U is). The 
ground state and the third excited state are widely separated, the interlevel spacing being 
proportional to 2^Jn for small E/T and U/T. 



A. Tripartite entanglement between the AQDs and the pool 

We will now turn to the issue of entanglement generation in the tripartite system. The 
tripartite system is described by the pure state 

1^) = ci\n + 1, 00) + C2|n, 01) + c^\n, 10) + Ci\n - 1, 11) (10) 
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where Cj(t) are complex coefficient satisfying ^^^^^ = 1- Although an entanglement 

„.ea.u.e i= not u„a,nbiguously defined tea tripartite hybrid .y.te„t (one ,utnt and two 
qubits) one can calculate a concurrence [27| 



cm) 



\ 



J2^rpi (11) 



where pi is the density matrix of a subsystem i with the two other subsystems traced out. 
The system is unentangled if this concurrence is zero, and a maximum of the concurrence 
corresponds to the entanglement maximum. We obtain from our full density matrix p = 
the three reduced density matrices 

Pi = \ci\'^\n + l){n+l\ + {\c2\^ + \c3\'^)\n){n\ + \c4\'^\n-l){n-l\, (12) 

P2 = (icip+ic2r)io)(oi + (ic3r+ic4r)ii)(ii, (13) 

P3 = (|ClP + |C3P)|0)(0| + (|C2P + |C4P)|1)(1|, (14) 

where pi, p2 and are the reduced density matrices of the pool of bosons, the left AQD 
and the right AQD, respectively. The resulting concurrence ( ITT]) reads 



cm) 



\ 



3 - 3^^ - 2 (|C2|2|C3|2 + \c2\W + + + \c3\W) ■ 

i=l 

(15) 

In Fig. 3 (a) we display an example of the time- dependent concurrence C for f//T = 0.2. 
The initial condition for Fig. 3 and the rest of the plots in the paper is Ci(0) = 1, which 
corresponds to two initially empty quantum dots. As expected, the time dependence of the 
concurrence is more complicated than for the bipartite case. However, it is again possible 
to identify an entanglement period and the maximum concurrence over that period. Note 
that the maximum concurrence is larger than 1, which is not surprising since our system 
comprises a qutrit — a quantum particle with a larger entanglement capacitance than a qubit. 
In Fig. 3 (6) we plot the maximum concurrence versus interaction U/T for various n. We 
observe similar tendencies as in the previous bipartite case. The maximum concurrence is 
suppressed by increasing either the interaction or the number of bosons in the pool. 

Figure 4 shows the dependence of the entanglement period tent on the bosonic interaction; 
it is rather different from the bipartite case (see Fig. 2(6)). Importantly, for small U/T the 
entanglement period can reach values orders of magnitude larger than in the bipartite case, 
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t U/T 

FIG. 3: (a) Time-dependent concurrence C of the tripartite setup for E /T = 0.01, n = 10 and 
interaction U /T = 0.2. (6) The maximum vahie of concurrence Cmax of the tripartite system versus 
interaction {//T for n = 4 (black), n = 10 (red) and n = 30 (blue) bosons in the pool. 

so that the system remains entangled over extended periods of time. We also note that for 
sufficiently large n the dependence of tent on the interaction scales approximately as 1/f/, 
so that tentU/T ^ const (see the inset in Fig. 4). Note that because of the units of tent, the 
product tentU/T is T-independent. 

Next we determine the bipartite entanglement in the system when one constituent of the 
tripartite system is traced out. Entanglement measures for mixed bipartite systems are well 
defined. In particular, we use the entanglement of formation and the negativity to quantify 
the entanglement in the bipartite system. 
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U/T 



FIG. 4: Entanglement period tent versus interaction for n = 4 (black), n = 10 (red) and n = 30 
(blue) interacting bosons in the pool. The product t^ntU /T versus U/T is shown in the inset. 

B. Bipartite entanglement between the AQDs 

To start with, we trace out the bosonic pool and consider a mixed two qubit state de- 
scribed by the reduced density matrix 

/ |cip \ 

|C2P C2C% 
C3C*2 |c3p 
V |C4|V 



Pi 



(16) 



written in a standard two-qubit basis. Entanglement of formation Ef can be defined in this 



case as 



2l| 



Efipi) = f 



1 + v/l-C(pi)^ 



(17) 



with f{x) = —x\og2X — {1 — x) log2(l — x), and the concurrence 



C(pi) = max{0, Ai — A2 — A3 — A4}. 



(18) 
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FIG. 5: (a) Time dependence of the entanglement of formation Ef for U/T = 0.2 and n = 10. (6) 
The maximum values of the entanglement of formation Emax (dashed curves) and the negativity 
^max (solid curves) versus the interaction U/T for n = 4 (black), n = 10 (red) and n = 30 (blue) 
bosons in the pool. 

Here, A, are the square roots of the eigenvalues of pipi in descending order, where pi = 
{ay ® o-y)pl{ay ® (Tj^ ) , aud pi is a complex conjugate of pi. After some algebra we obtain 
Ai = 0, A2,3 = |ci||c4| and A4 = 2|c2||c3|, which have still to be ordered. This allows us 
to calculate the concurrence and the corresponding entanglement of formation according 

to (dZD. 



An alternative measure of entanglement is the negativity 22 1 



^^{p^) = M]\^ . (19) 

Here, | |pi | |i is the trace norm of the partial transpose p\ of the bipartite mixed state. Essen- 
tially Af{pi) is just the absolute value of the sum of negative eigenvalues of . After some 
algebra we find that only the eigenvalue A = (|cip + |c4p — a/(|ciP — IC4P) + 4|c2p|c3p)/2 
of the partial transpose of the density matrix f|T6|l can be negative. 

The periodic time dependence of the entanglement of formation is presented in Fig. 5. 
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We note that the entanglement of formation and negativity have the same time dependence, 
however, they are different in amphtude, which is not surprising since the entanglement of 



formation and negativity are different entanglement measures [28|, |29|. We therefore analyze 
the dependence of the maximal entanglement of formation E^ax and maximal negativity 
Nmax on the interaction energy U /T. The results for different boson numbers n are shown 
in Fig. 5(6). 

An unexpected result is that the curves Emax{U) and Nmax{U) are nonmonotonic. Each 
curve possesses a maximum at a characteristic value oi U = Uc- The value of Uc depends 
on the number of bosons ra, more precisely, Uc is increasing with decreasing n; on the other 
hand, the Emax{Uc) and Nmax{Uc) are constant. The nonmonotonic behavior of the curves 
is probably related to the fact that for U /T the main entanglement capacity resides in 
the bosonic reservoir, which is traced out in the case at hand. Indeed, when just one ADQ 
is traced out instead, we observe a monotonic dependence (see next section). 



C. Bipartite entanglement between a single AQD and the pool 



Finally, we trace out one of the qubits and discuss entanglement of the hybrid qubit-qutrit 
system described by the reduced density matrix 



P2 





\r |2 
|Cl| 


C1C2 





\ 




C1C2 


1^2 1 


















C3C4 


\ 








C3C4 


|C4|V 



(20) 



In order to determine the negativity we again need to find negative eigenvalues of the partial 
transpose of the density matrix (!20l) . It turns out that only two eigenvalues can be negative, 
(|c2p - V|c2|^ + 4|c3|2|c4|2)/2 and A2 



namely Ai 



(|c3p- V|c3r + 4|ci|2|c2|2)/2. We 

^"1 (see 



can now calculate the negativity as an absolute sum of negative eigenvalues of \\p2 
previous section). 

The entanglement of formation can be also calculated by doing a convex expansion of the 
density matrix p2 |30j : 



(21) 
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whereby the pure states |\1/q,) and are 



I*,) = ^(ci|n + l,0) + C2|n,l)), (22) 



a 

= (C3|n, 0) + C4|n - 1, 1)) . (23) 

Here, a = |cip + |c2p and (3 = |c3p + |c4p, so that a + /3 = 1. One should note that the 
decomposition (l2Ti) is unique because of the prohibition against superpositions of states with 



different particle numbers. The resultant entanglement of formation is 30 1 



E = aE^ + PEf,, (24) 

where Ea(/3) is the entanglement of formation of the state |\1/q,) (|^'/3)) calculated with cor- 
responding concurrencies Ca — 2|ciC2|/a and C/? = 2|c3C4|//3. 

Figure 6 shows the entanglement of formation Ej versus time for a fixed interaction 
strength; the dependence of the negativity is similar up to a scaling factor. In Fig. 6(6) 
the maximum entanglement of formation E^ax and maximum negativity Nmax are plotted 
versus interaction U/T for different bosonic occupations. The behavior of these quantities 
is monotonic and in line with the U /T-dependence of the maximum concurrence shown in 
Fig. 3(6): Smaller boson numbers n in the pool favor larger values of entanglement. 



V. CONCLUSIONS AND DISCUSSIONS 

We have analyzed qubit-qubit, qubit-qutrit and qubit-qutrit-qubit systems which can 
be realized in a cold atomic set-up. Qutrits have bigger entanglement capacity and offer 
potentially more advantages for quantum computation than qubits. However, they are 
difficult to handle, as they involve more degrees of freedom. Qubit-qutrit entanglement has 



been recently observed for the first time in a photonic system [3l|]. Cold atomic qutrits, 
which emerge due to the particle conservation, are characterized by the number of bosons 
and their interactions; both parameters can be easily tuned experimentally. 

We have demonstrated that both parameters play a crucial role in generating the maxi- 
mal possible entanglement between the subsystems as well as for the longest entanglement 
period. For example, smaller number of particles and smaller interactions in general favor 
the generation of entanglement. Note that increasing n (n 1) leads to the disappearance 
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FIG. 6: (a) Time dependence of the entanglement of formation Ej for U/T = 0.2 and n = 10. (6) 
The maximum values of the entanglement of formation Emax (dashed curves) and the negativity 
^max (solid curves) versus interaction U/T for n = 4 (black), n = 10 (red) and n = 30 (blue) 
bosons in the pool. 



of entanglement as expected 23|. Systems with large interactions {U/T ^ 1) do not exhibit 
entanglement. Interestingly, the entanglement period of the tripartite system is orders of 
magnitude larger than that of the bipartite system. We have also found a simple relation 
between entanglement period and interaction for large number of bosons (n ^ 1). 



32| . i.e., entanglement 



Note that in this work we dealt with so-called mode entanglement 
of modes in the Fock representation. It was debated whether or not mode entanglement can 
be effectively used for quantum computatiom because superselection rules prohibit coherent 
superpositions of different particle number, which introduces severe limitations. However, it 
has been shown recently that there are ways to overcome this drawback, for example, mode 



entanglement can be used for dense coding 33 |. 



Our results may also shed lig ht on the role played by interactions in the problem of 



coupled harmonic oscillators 3J]. One can assume from our findings that entanglement 
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between oscillators strongly depends on the ratio of interaction versus coupling. Another 
interesting extension our results could be the study of entanglement between atomic Bose- 
Einstein condensates and immersed impurities 35, 3^. 
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